Caustics due to Negative Refractive Index in Circular Graphene p-n Junctions 
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We show that the wavefunctions form caustics in circular graphene p-n junctions which in the 
framework of geometrical optics can be interpreted with negative refractive index. 
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The possibility of negative refractive index in nature 
was first analyzed theoretically by Veselago[H which was 
followed by other important works Q. Evidence for such 
left-handed metamaterials has been demonstrated by mi- 
crowave experiments Q. 

Another prominent candidate for such materials might 
be the graphene as it was proposed recently by Cheianov 
et al. They studied the transmission of electrons 

through a plane p-n junction of graphene and showed 
that the optics of electron flow, in the framework of ge- 
ometrical optics, can be described by optical refraction 
with negative refractive index. The graphene itself pro- 
vides many other peculiar electronic properties owing to 
the close similarity between the dispersion relation of 
two-dimensional massless Dirac fermions and the low en- 
ergy electronic spectrum of graphene (for review see [f| ) . 
However, an easy way of tuning the refractive index by 
gate potential may open up further research directions 
in graphene physics. For example, graphene might be 
utilized to fabricate properly designed electronic lens. 

One of the interesting subjects in geometrical optics 
is the caustics. A caustic is an envelope of a family of 
rays at which the density of rays is singular. The caustics 
have been extensively studied in thepast (for review see, 
eg, the work by Berry and Upstill @)- As in Ref. [1] for 
plane p-n junction, we show theoretically that in circular 
p-n junction of graphene the caustics can also arise in 
the wavefunction pattern of electrons, and the curves of 
caustics can be calculated using the well-known Snell's 
law with negative refractive index. 

To this end, we consider the scattering of ballistic in- 
cident electrons in graphene shown in Fig. [T] To demon- 
strate the formation of caustics in circular p-n junctions, 
we used a simple gate potential V(r) = VoQ(R-r), where 
O denotes the Heaviside function. The assumption of 
such a sharp potential is valid provided that Af > d, 
where Xp is the Fermi wavelength outside the junction 
and d is the characteristic length scale in which the scat- 
tering potential varies. Moreover, to prevent interval- 
ley scattering it is neccessary that d a, where a is 
the lattice constant of graphene 0, 0] • The negative re- 
fractive index is realized by p-n junctions of graphene 
in which the energy E of the incoming electrons at the n 
side is positive (belonging to the conduction band), while 
at the p side of the junction the potential Vq is chosen 



FIG. 1: Incident plane wave of electron in single layer 
graphene is scattered by a rotational symmetric potential 
V(r) (circular region). 



as E < Vq, therefore the electrons belong to the valence 
band. Here we focus on this case (0 < E < Vq) but its 
generalization to arbitrary values of E is straightforward. 
Note that various properties of elastic electron scattering 
in graphene has already been studied [8(. Here we in- 
vestigate the wavefunction pattern inside the junction. 
Pronounced formation of the caustics can only be seen 
in the quasi-classical limit, ie, in case of R 3> A^, but no 
such condition is assumed in our exact calculation of the 
wavefunctions. 

The Hamiltonian for the above described scattering 
problem for energies E close to the Dirac point may be 
written as 



H = H + V(r)t = c er • p + V(r)l. 



(1) 



Here c is the Fermi velocity, p = —ihd/dr, while er = 
(a x , a y ) and 1 are the Pauli matrices and the unit matrix 
acting in isospin space. 

Consider the elastic scattering of incoming electrons 
governed by the Hamiltonian H . The scattering of the 
incident electron can be calculated by the generalization 
of the well-known partial wave method (see eg, ) . The 
wavefunction (in polar coordinates) describing the scat- 
tering of a single incoming partial wave h\ (defined be- 
low) outside the scattering region (r > R) can be writ- 
ten in terms of the cylindrical wave eigenfunctions of the 
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Hamiltonian H with energy E > 0: 
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while inside the scattering region (r < R) the wavefunc- 
tion reads 



(in) 



AjXj, where 
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Here the pseudo angular momentum j G J = 
{...,-§, §,...}, (hf>) is an outgoing (in- 

coming) cylindrical wavefunction corresponding to d = 1 
(d = 2). The wavenumbers are fc out = E/(hc) > and 
fc in = \E—Vo\/(he) > 0, while J n refers to the Bessel func- 
tion of the first kind, and H„' and are the Hankel 
functions of the first and second kind, respectively [loj . 
To construct the eigenfunctions we used the fact that for 
rotational symmetric potential [J Z ,H] = holds, where 
J z = —ihd v + Ha z /2 is the pseudo angular momentum 
operator. Therefore, the pseudo angular momentum is a 
conserved quantity in the scattering process. 

The scattering matrix Sj and the amplitude Aj are 
determined from the boundary conditions, ie, from the 
continuity of the total wavefunction at the boundary of 



the junction: \E^ out - > (r 
it is easy to find that 

Sj = 



R,(p) = lf n) (r = R,ip). Then, 
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where X = k out R and X' = k- m R. 

We now consider the scattering of an incident plane 
wave of electron in graphene for r > R. Such an eigen- 
state with energy E has the form 

O ik out r cos(ip—ipi) 



®<Pi(r, tp) 
vi<Pi) 




where (5a) 
(5b) 



and ifi denotes the direction of the propagation of the in- 
cident electron. Using the properties of the Hankel func- 
tions [Io| one can show that the partial wave expansion 



of $„ 
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Without the loss of generalization, we choose the direc- 
tion of propagation to be parallel with the x axis, which 
means ifi — in (|5a[) . Then the wavefunction for r > R 
is given by 



(out) _ 
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and the wavefunction for r < R is given by 

* (in3 =^E^'"H^- (7b) 

Note that the second term in (fTaj) is the scattered wave 
due to the scattering of the incident plane wave <f> Vi on 
the scattering region described by the potential V(r) — 
VoQ(R-r). The scattering cross section can be obtained 
from the asymptotic form (r 3> R) of the scattered wave. 

Equations ((3]) and |(7J) allows us to calculate exactly 
the wavefunctions both inside and outside the junction. 
Note that the wavefunctions depend only on the two di- 
mensionless parameters, fci n i? and k out R- Figures and [3] 
show how the incident plane wave (from direction ifi = 0) 
penetrates into the circular region of the p-n junction. 




FIG. 2: |>3'| 2 (in scale of logarithmic to base 10) is plotted 
inside and outside, close to the scattering region (the dashed 
line shows the boundary of the p-n junction). Here k[ n R = 300 
and koutR = 300 corresponding ton = — 1, and x and y are in 
units of R. The solid (dotted) line corresponds to the caustic 
for p — 1 (p = 2) (see the text). 

For r < R, the formation of the caustics in the wavefunc- 
tion patterns is clearly visible along the solid and dotted 
lines. Vodo et al. found a similar wave pattern exper- 
imentally by focusing a plane microwave with a plano- 
concave lens fabricated from a photonic crystal having 
negative refractive index [llj . 

In a sharp p-n junction of graphene the optics of elec- 
tron flow (4 1 is very much the same as in photonic crys- 
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FIG. 3: The same as in Fig. rjwith fe n i? = 300 and k out R 
200 corresponding to n = —1.5. 



tals with negative refractive index. In a similar way as in 
Ref. [4| for our circular p-n junction the Snell's law reads 



sin a 

sin/3 



kout 



(8) 




FIG. 4: The basic geometry of ray paths inside the circular p- 
n junction. The incident ray from left with impact parameter 
b and angle of incidence a is refracted at the boundary of the 
junction with angle f3 and then after p — 1 internal reflections 
exits from the junction. Here p — 3 and we set n — —1.3 and 
a = 60°. 




where a and /3 are the angle of incidence and refraction, 
respectively. Since in our calculation both fc; n and fc ou t 
are positive the refractive index becomes negative simi- 
larly as for properly designed photonic crystals. 

In what follows we show that the intensity maximum in 
the wavefunction patterns is around the caustic and can 
be understood by the Snell's law ([5]) with negative refrac- 
tive index. Figure 2] shows how the incident ray refracts 
at the boundary of the circular p-n junction and then af- 
ter p — 1 internal reflections exits from the junction. One 
can classify the different ray paths by the impact param- 
eter b = i?sinev and the number of chords p inside the 
circle corresponding to p— 1 internal reflections. Incident 
rays with varying impact parameters (— R < b < R) form 
a family of rays inside the circle. The envelope of this ray 
family results in a caustic as shown in Fig. [Sj Each of 
the chords has its own caustic. 

Using differential geometry to calculate the envelope of 
family of curves we find that the curve r c (in Descartes 
coordinates shown in Fig. rjj of the caustic of the pth 



FIG. 5: Caustic inside the p-n junction formed from the enve- 
lope of the refracted ray paths of different incident rays. The 
thick solid line is the caustic curve calculated from (JSJ. Here 
p — 1 and n = —1. 



chord is given by 
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Here a varies between —tt/2 and 7r/2 and the prime de- 
notes the derivation with respect to a. This is a para- 
metric curve for caustics with parameter a. Figures [5] 
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and [3] show the caustics calculated from Eq. ([9]). In both 
cases, one can see that the location of the caustics formed 
from the interference pattern of the exact wavefunctions 
inside the junction agrees very well with that obtained 
from SnelPs law with negative refractive index. The caus- 
tics for p > 2 are less visible since after each internal 
reflection the intensity of the rays decreases. 

The caustics in a circular p-n junction belong to the 
class of cusp according to the catastrophe optics [6(] . The 
location of the cusp (r,ip) — (r cusp ,7r) (in polar coor- 
dinates) for the pth chord can be obtained from ^ by 
setting a = and we find 



' cusp 



(p) 



(-If 



l + 2p 



R. 



(10) 



As can be seen in Fig. [5] the paraxial (a <C 1) incident 
rays entering into the circular region of the graphene 
junction are focused at the focal point for p = 1. The 
focal length / measured from the point (r, ip) = (R, ir) 
is given by the location of the cusp: f = R — |r cusp | = 
i?|n|/(|n| + 1). The same is obtained when the refractive 
index is replaced by its negative value in the expression 
of the image focal length defined in ordinary geometrical 
optics jl2| . To demonstrate the negative refractive index 
in the photonic crystal experiment ll| the focal length 
was measured. 

To realize the predicted caustic formation and focus- 
ing effect in an experiment, one needs to collimate a 
monodirectional electron beam onto the circular scatter- 
ing region. This might be achieved utilizing a smooth pla- 
nar p-n junction which is known to transmit only those 
quasiparticles that approach it almost perpendicular to 
the p-n interface 13[. Therefore, a possible experimental 
setup could be built up from a source electrode, a selec- 
tively transmitting smooth planar p-n junction, the cir- 
cular scattering region and a drain electrode. Under bias, 
the spatial dependence of the charge density of the trans- 
ported electrons might be measured by scanning probe 
techniques. 

Circular p-n junctions together with planar ones stud- 
ied earlier [4| can be a building block of electron optics 
in graphene. However, the refractive index varies with 
energy of the incoming electrons therefore the tempera- 
ture needs to be low enough for sharp images. As it is 
mentioned in the introduction, one needs experimentally 
that a -C d -C Aj? = 27r/fc out (for sharp potential barrier 
and absence of intervalley scattering). Another condition 
required to observe sharp interference patterns around 
the caustics is k m R 1 corresponding to the quasi- 
classical limit. For example, setting E = 40 meV, and 
Vo = 80 meV (using gate potential), R — 800 nm, all of 
these conditions are fullfiled since d ~ 10 nm 14J, imply- 
ing that the refractive index is n = — 1, while k m R = 50, 
and Xf — 10d. 



In summary, we calculated inside and outside a circu- 
lar p-n junction of graphene the scattered wavefunction 
of an incoming plane wave of electrons due to a circular 
symmetric step-like potential. We showed that the scat- 
tered wavefunction inside the junction has an interference 
pattern with high intensity maximum located around the 
caustics calculated from Snell's law with negative refrac- 
tive index. 
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